DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information
  1. Map-based cosmology inference with lognormal cosmic shear maps

    ABSTRACT Most cosmic shear analyses to date have relied on summary statistics (e.g. ξ+ and ξ−). These types of analyses are necessarily suboptimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative tomore » single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps – and to the extent that the convergence field can be described by the lognormal model – our map posteriors allow us to reconstruct all summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the σ8–Ωm plane by $$\approx\!{30}{{\ \rm per\ cent}}$$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $$S_8=\sigma _8\Omega _{\rm m}^{1/2}$$ directions, meaning map-based inference fails to significantly improve constraints on S8.« less
  2. Approximate Bayesian Computation applied to the Diffuse Gamma-Ray Sky

    ABSTRACT Many sources contribute to the diffuse gamma-ray background (DGRB), including star forming galaxies, active galactic nuclei, and cosmic ray interactions in the Milky Way. Exotic sources, such as dark matter annihilation, may also make some contribution. The photon counts-in-pixels distribution is a powerful tool for analysing the DGRB and determining the relative contributions of different sources. However, including photon energy information in a likelihood analysis of the counts-in-pixels distribution quickly becomes computationally intractable as the number of source types and energy bins increase. Here, we apply the likelihood-free method of approximate Bayesian computation (ABC) to the problem. We considermore » a mock analysis that includes contributions from dark matter annihilation in Galactic subhaloes as well as astrophysical backgrounds. We show that our results using ABC are consistent with the exact likelihood when energy information is discarded, and that significantly tighter parameter constraints can be obtained with ABC when energy information is included. ABC presents a powerful tool for analysing the DGRB and understanding its varied origins.« less
  3. Numerical relativity simulations of the neutron star merger GW190425: microphysics and mass ratio effects

    ABSTRACT GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $$\lesssim 0.083\mathrm{\, M_\odot }c^2$$ with peak luminosity of 1.1–$$2.4\times ~10^{58}/(1+q)^2\, {\rm {ergmore » \, s^{-1}}}$$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$$0.1 \mathrm{\, M_\odot }$$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.« less
  4. Observational window effects on multi-object reverberation mapping

    ABSTRACT Contemporary reverberation mapping campaigns are employing wide-area photometric data and high-multiplex spectroscopy to efficiently monitor hundreds of active galactic nuclei (AGNs). However, the interaction of the window function(s) imposed by the observation cadence with the reverberation lag and AGN variability time-scales (intrinsic to each source over a range of luminosities) impact our ability to recover these fundamental physical properties. Time dilation effects due to the sample source redshift distribution introduce added complexity. We present comprehensive analysis of the implications of observational cadence, seasonal gaps, and campaign baseline duration (i.e. the survey window function) for reverberation lag recovery. We findmore » that the presence of a significant seasonal gap dominates the efficacy of any given campaign strategy for lag recovery across the parameter space, particularly for those sources with observed-frame lags above 100 d. Using the Australian Dark Energy Survey as a baseline, we consider the implications of this analysis for the 4MOST/Time-Domain Extragalactic Survey campaign providing concurrent follow-up of the Legacy Survey of Space and Time deep-drilling fields, as well as upcoming programmes. We conclude that the success of such surveys will be critically limited by the seasonal visibility of some potential field choices, but show significant improvement from extending the baseline. Optimizing the sample selection to fit the window function will improve survey efficacy.« less
  5. AutoEnRichness: A hybrid empirical and analytical approach for estimating the richness of galaxy clusters

    ABSTRACT We introduce AutoEnRichness, a hybrid approach that combines empirical and analytical strategies to determine the richness of galaxy clusters (in the redshift range of 0.1 ≤ z ≤ 0.35) using photometry data from the Sloan Digital Sky Survey Data Release 16, where cluster richness can be used as a proxy for cluster mass. In order to reliably estimate cluster richness, it is vital that the background subtraction is as accurate as possible when distinguishing cluster and field galaxies to mitigate severe contamination. AutoEnRichness is comprised of a multistage machine learning algorithm that performs background subtraction of interloping field galaxiesmore » along the cluster line of sight and a conventional luminosity distribution fitting approach that estimates cluster richness based only on the number of galaxies within a magnitude range and search area. In this proof-of-concept study, we obtain a balanced accuracy of 83.20 per cent when distinguishing between cluster and field galaxies as well as a median absolute percentage error of 33.50 per cent between our estimated cluster richnesses and known cluster richnesses within r200. In the future, we aim for AutoEnRichness to be applied on upcoming large-scale optical surveys, such as the Legacy Survey of Space and Time and Euclid, to estimate the richness of a large sample of galaxy groups and clusters from across the halo mass function. This would advance our overall understanding of galaxy evolution within overdense environments as well as enable cosmological parameters to be further constrained.« less
  6. TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

    The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellarmore » kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within 1.2σ, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a ~4σ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within ~0.6σ<, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models.« less
  7. FLASH pilot survey: detections of associated 21 cm H  i absorption in GAMA galaxies at 0.42 < z < 1.00

    ABSTRACT We present the results of a search for associated 21 cm H i absorption at redshift 0.42 < z < 1.00 in radio-loud galaxies from three Galaxy And Mass Assembly (GAMA) survey fields. These observations were carried out as part of a pilot survey for the Australian SKA Pathfinder (ASKAP) First Large Absorption Survey in H i (FLASH). From a sample of 326 radio sources with 855.5 MHz peak flux density above 10 mJy, we detected two associated H i absorption systems, in SDSS J090331+010847 at z = 0.522 and SDSS J113622+004852 at z = 0.563. Both galaxies are massive (stellar mass $$\gt 10^{11}\, \mathrm{M}_\odot$$) and have optical spectramore » characteristic of luminous red galaxies, though spectral energy distribution fitting implies that SDSS J113622+004852 contains a dust-obscured starburst with star formation rate ∼69 M⊙ yr−1. The H i absorption lines have a high optical depth, with τpk of 1.77 ± 0.16 for SDSS J090331+010847 (the highest value for any z > 0.1 associated system found to date) and 0.14 ± 0.01 for SDSS J113622+004852. In the redshift range probed by our ASKAP observations, the detection rate for associated H i absorption lines (with τpk > 0.1 and at least 3σ significance) is $$2.9_{-2.6}^{+9.7}$$ per cent. Although the current sample is small, this rate is consistent with a trend seen in other studies for a lower detection rate of associated 21 cm H i absorption systems at higher redshift. We also searched for OH absorption lines at 0.67 < z < 1.34, but no detection was made in the 145 radio sources searched.« less
  8. Red Dragon: a redshift-evolving Gaussian mixture model for galaxies

    ABSTRACT Precision-era optical cluster cosmology calls for a precise definition of the red sequence (RS), consistent across redshift. To this end, we present the Red Dragon algorithm: an error-corrected multivariate Gaussian mixture model (GMM). Simultaneous use of multiple colours and smooth evolution of GMM parameters result in a continuous RS and blue cloud (BC) characterization across redshift, avoiding the discontinuities of red fraction inherent in swapping RS selection colours. Based on a mid-redshift spectroscopic sample of SDSS galaxies, an RS defined by Red Dragon selects quiescent galaxies (low specific star formation rate) with a balanced accuracy of over $$90{{\ \rmmore » per\ cent}}$$. This approach to galaxy population assignment gives more natural separations between RS and BC galaxies than hard cuts in colour–magnitude or colour–colour spaces. The Red Dragon algorithm is publicly available at bitbucket.org/wkblack/red-dragon-gamma/.« less
  9. WISE-PS1-STRM: neural network source classification and photometric redshifts for WISE×PS1

    ABSTRACT We cross-match between the WISE All-Sky and PS1 3π DR2 source catalogues. The resulting catalogue has 354 590 570 objects, significantly fewer than the parent PS1 catalogue, but its combination of optical and infrared colours facilitate both better source classification and photometric redshift estimation. We perform a neural network-based classification of the objects into galaxies, quasars, and stars, then run neural network-based photometric redshift estimation for the galaxies. The star sample purity and quasar sample completeness measures improve substantially, and the resulting photo-z’s are significantly more accurate in terms of statistical scatter and bias than those calculated from PS1 properties alone.more » The catalogue will be a basis for future large-scale structure studies, and will be made available as a high-level science product via the Mikulski Archive for Space Telescopes.« less
  10. The DESI N -body simulation project – I. Testing the robustness of simulations for the DESI dark time survey

    ABSTRACT Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogues to validate data analysis pipelines and identify potential systematics. We compare three N-body simulation codes, abacus, gadget-2, and swift, to investigate the regimes in which their results agree. We run N-body simulations at three different mass resolutions, 6.25 × 108, 2.11 × 109, and 5.00 × 109 h−1 M⊙, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the Dark Energy Spectroscopic Instrument (DESI) statistical error for $$s\ \gt\ 20\more » h^{-1}$$ Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of 2.1 × 109 h−1 M⊙ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $$20\ h^{-1}$$ Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1 per cent for k ≤ 10 h Mpc−1. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-N-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy–halo relation will be presented in a follow-up study.« less
...

Search for:
All Records
Subject
Modern experimental methods

Refine by:
Article Type
Availability
Journal
Creator / Author
Publication Date
Research Organization